A toy model of interlayer pair hopping 
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We extend the correlated electron model of Baskaran [Mod. Phys. Lett. B5, 
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\ 643(1991)] to the case of coupled layers. We show that the nature of the non-Fermi 
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liquid ground state leads to the absence of electron-like quasiparticles at the Fermi 
surface, thereby suppressing coherent transport of electrons between the coupled 
layers ("confinement"). On the other hand, motion of singlet pairs between the 
layers is not blocked. We also discuss how pair tunneling can be used to construct a 
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superconducting ground state. 



Typeset Using REVTEX 



1 



The study of correlated electronic systems has attracted considerable attention in recent 
times. In particular, the arguments advanced by Anderson and collaborators in favor of 
(a) failure of Fermi liquid theory in two dimensions (b) spin-charge separation ||, (c) 
irrelevance of single particle hopping between coupled Hubbard chains and/or planes || 
and (d) interlayer pair hopping || have evinced a lot of debate ]5J. Though most of these 
features are present in the one dimensional Hubbard model, the two dimensional model is 
yet to be solved. Therefore, it would be instructive to examine an interacting system of 
electrons which would exhibit these features in two and/or higher dimensions. Such a study 
would further our understanding of strongly correlated electronic systems. 

A step in this direction was taken by one of us || who proposed and solved a model of 
strongly interacting fermions in d-dimensions. The advantage of this model is that it can 
be solved exactly to show the breakdown of Fermi liquid theory and the emergence of spin- 
charge decoupling, thus illustrating two of the four features we listed above. In this letter, we 
extend this work to the case of coupled layers and examine the remaining two features, viz., 
single particle hopping between the coupled layers and interlayer pair hopping. We show 
that while the non-Fermi liquid ground state leads to a suppression of coherent one-electron 
motion between the coupled layers, motion of singlet pairs is not blocked. Both these results 
are obtained exactly. Finally, we discuss how long range superconducting order can be 
obtained in the presence of a residual intralayer pairing mechanism acting in conjunction 
with pair tunneling and how the scale of the superconducting transition temperature T c is 
governed by the interlayer hopping matrix element t±. 

We begin by writing the model hamiltonian of ref. ||, 

H = J2e k 4 CT c ka + J Y^S k .S k > , (1) 

fccr fcfc' 

where S k = c ka a a p c k p is the spin operator (S = |) in momentum space and e k is assumed 
h 2 k 2 

to be of the form — — for simplicity |7J . The primed sum indicates the absence of the k = k! 
term in the summand. Thus H a is a t — J model in momentum space. As shown in [6], H Q 
can be solved exactly in d-dimensions and the ground state \G) is given by 
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ig% = p s n 4, n 4 T 4i io> , (2) 

ki<\k\<k 2 |fe|<fci 

where P# is a singlet projection operator and {a k } is any ^-dependent spin configuration in 
momentum space in the singly occupied annulus of radii k\ and hi such that the condition 
Z)fci<|fc|<fc 2 51 = is satisfied. The radii k\ and ki are determined by the equations 

£k 2 — e fci — J 

E = E > 

fci<|fe|<fc F fe 2 >|A;|>A: F 

where fop is the Fermi wave vector of the noninteracting Fermi sea. The ground state is 
therefore a filled doubly occupied Fermi sea upto radius k\ ; between radii k\ and ki, it 
is singly occupied and the electrons inside this annulus form a total singlet in momentum 
space. Therefore, the momentum distribution function (at zero temperature) is given by 



1 if \k\ < ki 
\ if ki < \k\ < k 2 

As there is no discontinuity in the occupancy at kp, viz., Zk F = 0, the result implies a 
complete breakdown of Fermi liquid theory. (We note here that it is possible to restrict the 
range of the A;-space superexchange interaction to a shell of desired width about the Fermi 
surface, leading to an asymptotic failure of Fermi liquid theory as opposed to a complete 
breakdown as above.) Since the ground state is known exactly, this feature can also be seen 
by writing down the explicit form of the one-electron Green's function G(k,u). G(k,u) is 
defined in the usual manner as 

/+00 
e~ luJt ( Tc k} {t) 4 T (0) ) dt . 
-00 

We use the equation of motion method to evaluate G(k, to) exactly. The time dependence of 
Ck\ will be governed by commutators that occur in the equation of motion for G(k,u>) such 
as, 

3 J 

[ H, c feT ] = -(e k - — ) c fcT - J c kl S p ~ J c feT E S p ■ 

V V 



The expectation values of such commutators are evaluated with respect to one of the in- 
finitely many ground states. Since the ground state is a total singlet it follows that 

3 J 

(4 T [ H , °k\ ]> = ~(e* - -j) (4 T c fct> 

(Here, we note that it is not necessary for us to know the exact form of the ground state 
wave function. We only need to know that the ground state is a total singlet.) A similar 
line of reasoning will hold good for the higher order terms in the expression for c k ^{t). This 
simplifies our calculation of G(k,u) and we find that G(k,u) is given by 

G(k,u) = - t~ 3J-. 

iu - (e k - 

The pole that is present at a; = e k in the case of the non interacting system is absent and 
instead we now obtain a pole at u> = e k — The new pole has its residue contributed 
by one of the many (~ N) degenerate states. All these states are characterized by a double 
occupancy at point k and one unpaired spin elsewhere in the singly occupied /c-space annulus 
of radii k\ and k 2 . It is this which gives a potentially incoherent character to the spectral 
weight. For example, by choosing J in (1) to be dependent on momentum, it is possible 
to lift the degeneracy of these states and this would spread the pole to a truly incoherent 
background. This behaviour should be contrasted with the shifted pole structure in say, a 
mean field CDW or BCS state where the new poles essentially correspond to one or two new 
eigen states that are not incoherent. Thus it is clear that the shifting of the pole at k F is 
not a simple energy denominator effect or an energy shift effect. We can also see this by 
considering the motion of singlets, viz., the Green's function of a singlet pair S(k, k'; u). We 
define 

/+°° * ■ . 

(Tb kk ,(t)bl k ,(0)} dt , 

-oo 

where 

frL' = ( c 1t c l'i ~ c kl c k'l) ' 

is the creation operator of a singlet bond between k and k'. As in the case of G(k,u>), we 
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can evaluate S(k, k'; cu) exactly. We find that 

S(k,k';uj) = - j^— , 

iu - (e k + e k >) 

i.e., the pole gets shifted back and the two particle (singlet) Green's function has the same 
form as that of the non interacting system! This is essentially due to the singlet correlations 
induced by the /c-space superexchange term in H Q . Thus we see that while the pole at k F 
disappears in the case of the one-electron Green's function it persists in the two electron 
(singlet) Green's function. Obviously, this cannot happen if the vanishing of the pole were 
a simple energy denominator or an energy shift effect. 

We now ask if coherent hopping of electrons between coupled layers is possible in the 
presence of an interlayer hopping term 

£fj-(*0 ( C L dka + h.c. ) , 

krj 

Here, t±(k) is the dispersion along the c-axis and c and d represent the coupled layers. 
(Though our results are valid for d-dimensions, we shall continue to use the term "layer" 
throughout.) The hamiltonian H is now given by 

+ E k *t±(k) (4,4. + h.c.). 
We check for c-axis current by evaluating the Green's function 

G\k,u) = -i ( Tcl(t)d ka (0) )e~ wt dt . 

J —CO 

If this has a pole structure, it would mean that electrons can be transported along the c- 
direction. In such a case, the interlayer coupling would be termed a "relevant operator". 
We evaluate G ± (k, cu) by defining the operators a ka and (5 ka , 

Otko = (Cjk CT + d ka ) 

fiko = —7= (Cfco- — dka) , 
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and reexpressing the hamiltonian in terms of these operators. It is easily verified that H 
can be rewritten as 

+ 77 D fc - D fc' "+" D fc "+" •<-V 



J 



where we have defined 



+ ? E (^ + ^ Q )-(^ + ^ Q ) ; 
z kk' 



1 t 

'S'fc = 2 <V afc ^ ' 



qa(3 



1 



77 «L ^ Pkv , etc. 



2 J fe^t u nv 

Since ( Tc|. (T (t) c? fccr (0) ) = ( Td\. a (t) c ka (0) ) by symmetry, we can write 



r+°° 

J — oo 



^ — ILOt 



( Talit) a ka (0) > - < Tpl(t) (3 k<7 (0) ) 



Though we could not write down the exact ground state wave function we have nevertheless 
been able to obtain the exact form of G ± (k, uo) by exploiting the symmetry between the two 
layers c and d and by using the fact that the ground state wave function \GS) must be a 
total singlet, viz., 

£(ST + S*") \GS) = 0, 
k 

(where S c ~ and S d ~ are the usual spin lowering operators) as well as the symmetry between 
the coupled layers by virtue of which 

J2S c k z -S d k z \GS) =0 etc., 

k 

In terms of the a and the (3 operators, we then have the following conditions satisfied by 
the ground state \GS). 

E( s r+ 5 ?") \gs) = o, 

k 

E (s** + st) \GS) = 

k 
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and 

EiSf' + Stl \GS) = 0, 

k 

E (sf z + si az ) \GS) = o. 

k 

Using these conditions, we arrive at the form of the Green's function for interlayer motion, 

„t ™. \ I R. flt 



G (k,u) oc 



("L a to) (PkaPha) 



iw - (e k - t ± (k) - f ) zu; - (e fc + - f ) 

Since in general, (a^o;^) and (fl ka (3 ka ) < 1, we see again that there will be no pole structure 
in the Green's function for energies less than ^ as in the single layer case. As discussed 
before, these two poles are also potentially incoherent and modifying the form of J would 
spread the pole to an incoherent background. Thus the non Fermi liquid nature of the 
ground state persists in the presence of interlayer hopping and this causes the bare hopping 
matrix element to get renormalized to zero. This is the phenomenon of "confinement" which 
is being studied extensively in coupled Hubbard chains/planes. 

We now turn our attention to the Green's function for singlet transport between the 
layers S cd (k, k'; u). It is not difficult to see that S cd (k, k'; u) is given by 



S (k, k'\ uj) oc 



iuj - (e k - t±(k) + e k > - t ± (k')) iu - (e k + t±(k) + e v + t ± (k')) 

The pole structure is got back when singlet transport is considered. So we see that while one 
electron motion between the coupled layers is blocked (the corresponding Green's function 
being completely incoherent), singlets can be transported freely. Another way of interpreting 
this result is to say that while the state describing an electron-hole excitation between the 
layers, c\, a d ka \N c ; N d ) suffers an orthogonality catastrophe 

((N + l) c ;(N-l) d I c{ a d ka \N c ;N d ) = 0, 

there is a finite overlap between states that are connected by pair hopping, viz., 

((N — 2) c ; (N + 2) d | b d k , b c kk , \N C ; N d ) £ 0. 



Now that we have shown that it is possible to transport pairs of electrons between the 
two coupled layers, the next obvious question would be to ask if such a process leads to 
superconductivity. Since we have calculated the interlayer two particle Green's function, 
we can, in principle calculate the Josephson current between the coupled layers and discuss 
the nature of the superconducting state. But we adopt the following perturbative approach 
which not only illustrates the pair hopping mechanism but also enables us to compare our 
results with those obtained earlier in the context of high temperature superconductivity from 
interlayer tunneling |§], ||, ||10|| . As we have shown that the interlayer single electron Green's 



function does not have a pole at whereas the two electron (singlet) Green's function does 
, it is enough to consider only pair hopping processes that are generated perturbatively. Let 
us consider the interlayer hopping as perturbing the initial ground state \G) Q by breaking, 
for instance, a singlet bond between momentum states k and —k with the rest of the wave 
function being undisturbed. This results in the creation of a particle hole excitation between 
the layers which we treat as an intermediate state in perturbation theory. It is easy to see 
that second order perturbation theory leads to an effective interaction of the form 

- 37 X^J-( fc ) c h c -ki d ~kid k i ■ (3) 

4 k 

This hamiltonian resembles the holon pair hopping hamiltonian that was first introduced 
by Wheatley, Hsu and Anderson j|. However there is an important difference since the 
pair hopping interaction above is extremely local in fc-space. This is because the interlayer 
hopping matrix element t±(k) (which generates the pair hopping process) conserves the 
momentum parallel to the layers. The significance of this fact and its implications were first 
discussed by Anderson M. Such a pair hopping hamiltonian has been studied recently by 



several people ||, JT0|. In particular, one can for instance do a mean field calculation with the 
pairing interaction as given by (3) and obtain the superconducting transition temperature 

r m t±{k)max 
kB+c ~ J • 

Note that T c depends on the pairing interaction in a simple manner unlike the BCS case. 



But the extreme (fc-space) locality of the pairing interaction will suppress any finite temper- 
ature phase transition since the phases of various pairing amplitudes (cLcl fc |) do not get 
correlated. However, it is clear that any residual pairing interaction such as an intralayer 
BCS interaction which is non local in fc-space will correlate the phases of various pairing 
amplitudes and will result in long range superconducting order. It should be emphasized 
that the scale of T c will still be set by the strength of the pair tunneling interaction rather 
than that of the intralayer interaction as long as -y > Vbcs I PL Q|| - 

Before we conclude we would like to point out that the model under consideration was 
introduced to mimic the single occupancy constraint in a realistic model such as the large 
U Hubbard model, where it is speculated that forward scattering of a pair of electrons 
becomes singular asymptotically close to the Fermi surface. The model discussed in this 
letter does not show this asymptotic behaviour but leads to a complete breakdown of Fermi 
liquid theory. We believe this model can be made more realistic by introducing a suitable 
momentum dependent J in (1). This is a topic of future study. 

To summarize, we have been able to illustrate several features of interlayer pair hopping 
that have been proposed in connection with high temperature superconductivity such as 
"confinement" , fc-space locality of the pair hopping interaction and the analytic dependence 
of T c on the strength of the pair hopping interaction. Our results assume added significance 
owing to the fact we have derived these features from a microscopic hamiltonian, albeit a 
toy one. 
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